Theory of superconducting fluctuations in the thinest carbon nanotubes 
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The low-energy electronic Hamiltonian for the thinest zigzag carbon nanotube, embedded into a 
dielectric host, is derived and its phase diagram is discussed. The specific multi-band structure and 
the microscopic form of the electron-electron interaction in this systems is considered. The interband 
repulsive interaction, which is almost unscreened, leads to a retarded intraband attraction between 
the electrons and can stabilize the Cooper pairs. For a dielectric constant of the host ed ~ 2 — 4, the 
theory predicts that the superconducting fluctuations should develop, which is in agreement with 
the experiment. For > 8, the density wave fluctuations should be amplified. Between the two 
phases, there is a metallic state where all two-particle fluctuations are suppressed. 

PACS numbers: 71.10.Pm, 73.63.Fg, 74.20.Mn, 74.70.Wz 



Carbon nanotubes fl attract the attention because of 
their unique geometrical structure, mechanical, chemi- 
cal or electrical properties, and possible applications in 
electronic devices [pi 9], A single- wall carbon nanotube 
(SWNT) is made of a graphite layer rolled up into a 
cylinder || with the smallest diameter about AA jjj. 
In such thinest SWNTs, superconducting fluctuations 
(SCFs) were observed with the mean-field critical tem- 
perature about 15 K |^| . Since it is a one-dimensional 
(1-d) system, the true off-diagonal long-range order is 
suppressed. Nevertheless, the SCFs might develop and 
the question arises about their microscopic origin. 

The low-energy electronic excitations in the metallic 
SWNT with a moderate diameter are described by the 
Luttinger liquid theory with four collective modes ||. 
The SCFs induced by the backscattering of the electrons 
might appear only at the lowest energy scale and it is 
rather unlikely that they can be ascribed to the occur- 
rence of such fluctuations in the 4A-SWNTs f§ . Theories 
considering lattice vibrations require a large reduction of 
the the repulsive electron-electron interaction to stabilize 
the Cooper pairs . The renormalization group analysis 
confirms that the SCFs due to phonons in the SWNTs 
can easily be destroyed by the competing Coulomb repul- 
sion H . The theory |?J explains the superconductivity in 
ropes of the SWNTs ||, where also the intertube hop- 
ping plays the role. However, in the experiment S, the 
SWNTs were embedded into a zeolite dielectric host and 
they were isolated from each others. Therefore, it seems 
rather unlikely that the SCFs in such a system may de- 
velop due to phonons. 

In the present paper we derive a low-energy electronic 
Hamiltonian for the thinest zigzag SWNTs embedded 
into a dielectric host. Because of a large curvature of 
the cylindrical shell in such SWNTs, there are three 1-d 
bands crossing the Fermi level. We find a phase diagram 
determined by the slowest decaying in space two-particle 



correlator. For realistic parameters, corresponding to the 
experimental realization (H), the SCFs are developed in 
the system. They are induced by the unscreened inter- 
band Coulomb interaction. Changing the dielectric prop- 
erties of the host, other possible phases are predicted. 

The simplest tight-binding approximation predicts 
that the (n, m) SWNT is metallic if mod(n — m, 3) = 
(where (n, m) is a chiral vector describing uniquely the 
infinitely long SWNT) ||. In such a case there are two 
bands crossing the Fermi level. 

When the curvature of the cylinder is not neglected, 
it turns out that from all of the metallic SWNTs only 
those with the (n, n) chiral vectors (armchair SWNTs) 
remain metallic and the rest are semiconductors with a 
small gap in the spectrum JhJ. However, for the thinest 
SWNT (n, 0) (zigzag SWNTs) with n = 4, 5 and 6, it 
has been shown that the Koster-Slater parameterization 
in the tight binding approximation is not adequate be- 
cause it does not capture correctly the hybridization be- 
tween the o~ — 7r orbitals |ll| . The functional density cal- 
culations within the local-density approximation (LDA) 
have shown that such nanotubes are metallic. Because of 
the very strong a — n hybridization, the nondegenerate 
band originating from a T — M line in the hexagonal Bril- 
louin zone is strongly repulsed and pushed down crossing 
the two-fold degenerate valance bands. As a result the 
Fermi level is lowered and there are three 1-d bands that 
are partially populated by the electrons leading to the 
metallic band structure. 

The low-energy theory for the thinest zigzag SWNTs is 
given by the three-band Tomonag a- Luttinger model with 
one nondegenerate band and two-bands that are two- fold 
degenerate. The corresponding 1-d Hamiltonian has the 
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following kinetic part 

2 r 

H o = -'EB^H / dx i>i jtr (x)d x i> rj Ax), (1) 

j=0 rr,r J 

where ip r ja {%) are slowly varying Fermi fields for the right 
r = R and the left r — L moving electrons with the 
spin cr in one of the three bands labeled by j. The non- 
degenerate band is labeled by j = and has the Fermi 
velocity v F . The two- fold degenerate bands are labeled 
by j = 1 and 2 and have the same Fermi velocity vf- 
The Fermi velocities might be determined from the LDA 
band structures, which were calculated in the supercells 
geometry of the SWNTs In the experimental real- 
ization H , the SWNTs were in 1-d channels of a porous 
zeolite (AIPO4-5) and, in principle, interacted with the 
walls of the host by the van der Waals forces. The quali- 
tative structure of the bands should not be changed with 
respect to the LDA results, however, small differences in 
the position of the Fermi level and the bottom of the 
r — M band might be expected. Therefore, the numeri- 
cal values for the Fermi velocities require more detailed 
investigation considering the presence of the host. In the 
following we keep these two velocities as parameters in 
the model. For further convenience we introduce the ra- 
tio i] = v f /vf and parameterize the model by vp and 

7]. 

The long-range Coulomb interaction between the elec- 
trons is modified by the zeolite host, which has a dielec- 
tric constant e^. One-dimensional wire embedded in a 
host with a different dielectric constant was studied in 
|l2"| where the potential describing the electron-electron 
and the electron-image-charge interactions was derived. 
Using now the wave functions centered at the cylindrical 
shell of the SWNT, we calculate the Fourier component 
V(q) of the electrostatic potential in the tube embedded 
into the dielectric host. Explicitly, the potential takes 
the form 



V(q) = — [h{qR)K Q {qR)- 
irei 



1 - f*) Ki(qa)K (qa) 



^K^qajl^qa) + K (qa)h(qa) 



(2) 



where a is the radius of the cavity in the zeolite, R < a 
is the radius of the SWNT, ei is the dielectric constant 
of the SWNT, and /„ and K n are modified Bessel func- 
tions. For small q the potential diverges as ln|gi?|. For 
large q the tail of the potential V(q) is suppressed be- 
cause of the image charges. Therefore, a scattering with 
large momentum transfer, in particular the backscatter- 
ing with q — 2k f , in the SWNTs embedded into the di- 
electric host should be suppressed and change the phase 
diagram only at the lowest energy/temperature scales. 

The low-energy Hamiltonian requires the screened 
electron-electron Coulomb potential obtained by remov- 



ing self-consistently high-energy degrees of freedom. In 
the random-phase-approximation the particle-hole polar- 
ization diagrams are only kept and it leads to the dynam- 
ically screened electron-electron Coulomb potential given 
by the algebraic equation J2u> zil'nn' (q, u)Vfi, mm , (q, u) = 
Vnn' mm' (q), where V nn - mm > (VJf, mm ,) is the matrix el- 
ement for the bare (screened) Coulomb potential and 
e«W = <*in<W - TLll'Viunn' is the matrix dielec- 
tric function expressed by the particlc-holc polariza- 
tion function ILj'(g, u>) fl3[ ]. The indexes correspond 
to the Bloch functions of the corresponding band. In 
the static limit lu = and for the momentum trans- 
fer q — > we obtain the screened intraband poten- 
tial V^ooq « Voooo(l - 2Vi m riii )/e intra and V£ ln = 

V2222 ~ ^llll(l — MjOOOnpo)/ Eintra with Cintra ~ (1 — 

VboooIIooXl — 2ViinIIii) Q. However, the interband 
Coulomb potential remains almost unscreened within 
this limit, i.e V S 101 = V^ 202 ps V i i/e inter with e mter w 
1. 

The effective interaction between the electrons is de- 
scribed by the two-particle part of the Hamiltonian 



Hp = 



iEEE 

i,j=0 rr' 00' 



dxdx' p ri <y{x)U i: j{x - x')p r > ja >(x'), 

(3) 



where p r ia(x) — iplj a (x)ilj r j a (x) is the density fluctuation 
operator. Taking the bare potential Vuu(q) — Vijij(q) = 
V{q) and approximate it by its average value V av calcu- 
lated between the infrared-cutoff gj, nf = 2ir/L, where L 
is the nanotube length, and the ultraviolet-cutoff g" 1 * = 
5 • 10 8 to -1 (which gives the energy cutoff E c = 0.1 eV), 
we obtain the following local potential Uij(x — x') = 
UijS(x - x') with C/oo = Kv/[1 + 2V &v /TTr,v F ] = U , 
U n = U22 = U12 = Kv/[1 + War/imp] = Ux, and 

U01 = U02 - Kv = w m. 

To find exponents in the one- and two-particle correla- 
tion functions in the thinest SWNT, we need the explicit 
solution of the Hamiltonian H = H Q + Hp, where H 
and Hp are expressed by Eqs. (@) and (|). This can be 
achieved by using the nonabelian bosonization where con- 
jugate fields (j>j a (x) and Hj a (x) are introduced for each 
of the band j and the spin a Jlq] . The Hamiltonian sep- 
arates into different sectors if we introduce the collective 
variables 
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for nondegenerate and degenerate bands, respectively. 
There are six types of the collective oscillations in the 
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system: four of them are neutral and two of them are 
charged. This feature distinguishes the thinest zigzag 
SWNTs from the other metallic SWNTs with moderate 
diameters. 

The interaction appears only between the charged 
density modes. Rescaling the fields <f> = <p/vK and 
n = VKT1, we obtain the following Hamiltonian in the 
charge sector 



rir = I <l.r I -f- \ H,V:ri + d x cj> 0c (x) 



TABLE I: Exponents for the two-particle correlation func- 
tions. 



Tl 2 c (x)+ d x <j> c {x) 



+ A 



d x <j} 0c {x) d x 4> c (x) ,(6) 



where the Luttinger liquid parameters are Kq c 



Uo/nvp and K c 



lA/1 



renormalized velocities are v Fc 



\- 2Ui/ttvf, and the 
Vp/K c0 and v Fc = 
vf/K c . The coupling between cj) 0c and (j) c modes is given 
by A = U'^2K 0c K c /ir. The particular band structure 
of the thinest SWNTs leads to the existence of the two 
charge modes which are coupled by the Coulomb inter- 
action. 

The Hamiltonian (^|) is bilinear and can be diagonal- 
ized taking the free-boson form with two characteristic 
velocities «£. = [v%, 2 +v$±y/(v$' 2 - v 2 ) 2 + AX 2 d F i F }/2 
0. The solution is real only if v± c > 0, i.e. v f vf > A 2 . 
For large A the system is unstable toward the long-wave 
length density fluctuations. This kind of instability is 
known as the Bardeen-Wentzel instability (BWI), origi- 
nally discussed in the case of the 1-d electrons coupled 
to phonons |l9). The BWI for two electrostatically 
coupled Luttinger liquids was also discussed in Ref . [jl9| . 
In a proximity to the BWI but still on the metallic side, 
the repulsive interband interaction can lead to a retarded 
intraband attraction and the electrons may form Cooper 
pairs. This is a possible microscopic mechanism respon- 
sible for developing the SCFs in the 4A-SWNTs. 

In the neutral sector, the Hamiltonian takes the sim- 
ple free-boson form with four characteristic velocities 
v Fs = v f /Kq s and vfu = vf/K u , where Kq s = 1 and 
K v = 1, and v — s,f or x. We have left the Luttinger 
liquid parameters (Kq s and K v ) explicitly because they 
are expected to flow due to a renormalization caused by 
backscattering processes. 

The one-particle correlation functions at zero temper- 
ature decay with a distance (^ r j (T (x)^ - (T (0)) ~ x~ aj , 
where the interaction dependent exponents ctQ — K$ s /A-\- 
l/AK 0s + K 0c A /4 + B /4K 0c , and a x = a 2 = K s /8 + 
l/8K s +K f /8+l/8K f +K x /8+l/8K x +K c A/8+B/8K c . 
The coefficients As and Bs are given by certain combina- 
tion of the effective velocities and Bogoliubov parameters 

Deviations form ctj = 1 indicate the Luttinger liquid 
behavior and can be experimentally detected by mea- 
suring the tunneling conductance through the SWNT 
p?j|. Each transport channel (j) contributes in paral- 
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lei to the total conductance G by its partial conduc- 
tance Gj(T,V sd ) = T a '- 1 F{eV sd /kT), where T is the 
temperature, V s d is an applied bias voltage, and F(z) 
is a known universal function p0[ . We predict the fol- 
lowing form of the conductance G(T, V sd ) = (T Qo_1 + 
2T ai - 1 )F{eV sd /kT) in the thinest zigzag SWNTs. 

Two-particle correlation functions at zero temperature 
decay according to the power law x~^ j (with (3j = 2 in 
the noninteracting case). The slowest decaying correla- 
tion function indicates which kind of fluctuations dom- 
inates in the system. We calculate the exponents (3j 
corresponding to: a singlet superconducting order pa- 
rameter (SSC), a triplet superconducting order parame- 
ter (TSC), a charge/spin density wave order parameter 
(CDW/SDWz), and a spiral density wave order param- 
eter (SDWs), with the nesting vectors Qj = 2k F for the 
CDW/SDWz and SDWs cases (k F is the Fermi vector in 
the corresponding band). The results for the exponents 
of the intraband correlation functions are summarized in 
the Table I. The interband two-particle correlation func- 
tions decay much faster with x and are not relevant here. 

In Fig. [I] we present a phase diagram at zero tempera- 
ture in the parameter space defined by the average bare 
interaction V" av (in units o{Hvf) and the ratio r\ = v f /vf- 
For large 77 (> 6) and small V av /vF the density wave fluc- 
tuations (DWFs) in the j = 1 and 2 bands appear to have 
the longest decaying length. For moderate V av /hvF, the 
DWFs dominate in the j = band and finally there is a 
wide part of the diagram where the SCFs develop in the 
j = 1 and 2 bands. This phase is terminated by the line 
where the BWI occurs. For moderate 77, an additional 
phase appears where all kind of the two-particle fluctua- 
tions are suppressed (the gray area in Fig. ^) . For small 
77 and Vav/hvF, the DWFs evolve in the j = band and 
for larger V av /hvF, close to the BWI line, the SCFs are 
amplified in the j = band. 

The SCFs with the longest decaying length develop in 
a band with a smaller renormalized velocity. In other 
bands, with a larger renormalized velocity, the SCFs are 
also amplified (with (3j < 2). However, when the SCFs 
are enhanced, the DWFs are suppressed (with flj > 2) 
in all of the bands. Since K$ s = K v = 1 , the exponents 
in the SSC and the TSC correlation functions are the 
same. This is also the case for the CDW/SDWz and the 
SDWs correlation functions. Further analysis, includ- 



4 



3 - 



cdw/sdw 



e. = 40,20, 10, 




1.4 



unstable 



In conclusion, fabricating the SWNTs inside the zeolite 
hosts with different dielectric constants would provide a 
rare opportunity to control the electron-electron inter- 
action and to examine the different phases in such 1-d 
systems. 
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FIG. 1: The phase diagram for the thinest zigzag carbon nan- 
otubes in a parameter space given by the average bare interac- 
tion Vn V /h,VF and the ratio r\ = v%/vf of the Fermi velocities. 
The sector of the diagram with the dominant superconducting 
fluctuations in the j = 1 and 2 (j — 0) bands is marked by SC 
(SCO) and the portion with the dominant charge/spin density 
wave fluctuations in the j = 1 and 2 (j — 0) bands is indi- 
cated by CDW/SDW (CDW0/SDW0). The gray area shows 
the metallic phase where all of the two-particle fluctuations 
are suppressed. The thick solid line points out the Bardeen- 
Wentzel instability^ Large dots represent the averages for the 
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